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Abstract. We classify all pairs of recurrence relations in which two Hahn or dual Hahn 
polynomials with different parameters appear. Such couples are referred to as (dual) Hahn 
doubles. The idea and interest comes from an example appearing in a finite oscillator 
model [Jafarov E.I., Stoilova N.I., Van der Jeugt J., J. Phys. A: Math. Theor. 44 (2011), 
265203, 15 pages, arXiv:1101.5310]. Our classification shows there exist three dual Hahn 
doubles and four Hahn doubles. The same technique is then applied to Racah polynomials, 
yielding also four doubles. Each dual Hahn (Hahn, Racah) double gives rise to an explicit 
new set of symmetric orthogonal polynomials related to the Cliristoffel and Geronimus trans¬ 
formations. For each case, we also have an interesting class of two-diagonal matrices with 
closed form expressions for the eigenvalues. This extends the class of Sylvester-Kac matrices 
by remarkable new test matrices. We examine also the algebraic relations underlying the 
dual Hahn doubles, and discuss their usefulness for the construction of new finite oscillator 
models. 
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1 Introduction 

The tridiagonal (N + 1) x (N + 1) matrix of the following form 

\ 


Cn +i — 


1 


/ o 

N 0 


N - 1 0 


2 0 N 

1 0 / 


( 1 . 1 ) 


appears in the literature under several names: the Sylvester-Kac matrix, the Kac matrix, the 
Clement matrix, .... It was already considered by Sylvester [28], used by M. Kac in some of his 
seminal work [17], by Clement as a test matrix for eigenvalue computations [9], and continues 
to attract attention [6, 7, 29]. The main property of the matrix C/v+i is that its eigenvalues are 
given explicitly by 


-N, -N + 2, -N + 4,..., N -2, N. (1.2) 

Because of this simple property, Cjv+i is a standard test matrix for numerical eigenvalue com¬ 
putations, and part of some standard test matrix toolboxes (e.g., [12]). 




2 


R. Oste and J. Van der Jeugt 


One of the outcomes of the current paper implies that Cjv+i has appealing two-parameter 
extensions. For odd dimensions, let us consider the following tridiagonal matrix 


C , 2at+i(7, <5) 


/ 0 2 7 + 2 
21V 0 2 

28 + 2N 0 27 + 4 

2N -2 0 


4 


\ 


V 


25 + 4 0 27 + 21V 

2 0 2 N 

28 + 2 0 / 


(1.3) 


In the following, we shall sometimes use the term “two-diagonal” [10] for tridiagonal matrices 
with zero entries on the diagonal (not to be confused with a bidiagonal matrix, which has also 
two non-zero diagonals, but for a bidiagonal matrix the non-zero entries are on the main diagonal 
and either superdiagonal or the subdiagonal). So, just as C 2 N +1 the matrix (1.3) is two-diagonal, 
but the superdiagonal of C 2 N+ 1 , 


[lj 2, 3,4,..., 21V — 1, 2N] 


is replaced by 


[2 7 + 2, 2, 2 7 + 4,4,..., 2 7 +21V, 21V], 
and in the subdiagonal of C 2 N+ 1 , 

[21V, 21V- 1,21V-2,..., 3, 2,1] 
the odd entries are replaced, leading to 

[21V, 28 + 21V, 21V — 2,... ,28 + 4, 2, 2(5 + 2], 

Clearly, for 7 = 8 = — ^ the matrix C^tv+i^, <5) just reduces to C^jv+i- One of our results is 
that C , 2 jv+i( 7 , 8) has simple eigenvalues for general 7 and 8, given by 

0, ± 2 v / l ( 7 + <5 + 2 ), ±2^2(7 +-5 + 3), ±2^3(7 + (5 + 4),..., ±2-//V ( 7 + <5 + N + 1 ). 

This spectrum simplifies even further for 8 = —7 — 1; in this case one gets back the eigenva¬ 
lues ( 1 . 2 ). 

For even dimensions, we have a similar result. Let C 2 n{'Y, <5) be the (21V) x (21V) tridiagonal 
matrix with zero diagonal, with superdiagonal 

[2 7 + 2, 2, 2 7 + 4,4,..., 21V - 2, 2 7 + 21V] 

and with subdiagonal 

[2(5 + 21V, 21V - 2, 28 + 21V - 2,..., 4, 28 + 4, 2, 28 + 2], 

Then C 2 n(+, 8) has simple eigenvalues for general 7 and <5, given by 1 

± 2^(7 + 1)(<5 + 1), ±2^/(7 + 2)(<5 + 2),..., ±2^(7 + N)(S + N). (1.4) 


1 In this paper, we shall usually take 7 > — 1 and <5 > — 1, yielding real eigenvalues for the matrices under 
consideration. But the expressions for the eigenvalues remain valid even when these conditions are not satisfied. 
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This spectrum simplifies for S = 7 , and obviously for 7 = 5 = — \ one gets back the eigenva¬ 
lues ( 1 . 2 ) since in that case 62 ^( 7 , <5) just reduces to C 27 V. 

What is the context here for these new tridiagonal matrices with simple eigenvalue proper¬ 
ties? Well, remember that Cjv+i also appears as the simplest example of a family of Leonard 
pairs [24, 30]. In that context, this matrix is related to symmetric Krawtchouk polynomials [13, 
19, 23]. Indeed, let K n (x) = K n (x;^,N), where K n (x;p,N) are the Krawtchouk polyno¬ 
mials [13, 19, 23]. Then their recurrence relation [19, equation (9.11.3)] yields 

nK n -i(x) + (N — n)K n+ i(x) = (N — 2x)K n (x), n = 0,l,...,N. (1.5) 

Writing this down for x = 0,1,..., iV, and putting this in matrix form, shows indeed that the 
eigenvalues of Cn+i (or rather, of its transpose C^ +1 ) are indeed given by (1.2). Moreover, it 
shows that the components of the kill eigenvector of Cj^ +1 are given by K n (k). 

So we can identify the matrix Cn+i with the Jacobi matrix of symmetric Krawtchouk poly¬ 
nomials, one of the families of finite and discrete hypergeometric orthogonal polynomials. The 
other matrices CV( 7 , J) appearing in this introduction are not directly related to Jacobi matri¬ 
ces of a simple set of finite orthogonal polynomials. In this paper, however, we show how two 
sets of distinct dual Hahn polynomials [13, 19, 23] can be combined in an appropriate way such 
that the eigenvalues of matrices like CV(7, S) become apparent, and such that the eigenvector 
components are given in terms of these two dual Hahn polynomials. This process of combining 
two distinct sets is called “doubling”. We examine this not only for the case related to the 
matrix Cxi'}, $), but stronger: we classify all possible ways in which two sets of dual Hahn poly¬ 
nomials can be combined in order to yield a two-diagonal “Jacobi matrix”. It turns out that 
there are exactly three ways in which dual Hahn polynomials can be “doubled” (for a precise 
formulation, see later). By the doubling procedure, one automatically gets the eigenvalues (and 
eigenvectors) of the corresponding two-diagonal matrix in explicit form. 

This process of doubling and investigating the corresponding two-diagonal Jacobi matrix can 
be applied to other classes of orthogonal polynomials (with a finite and discrete support) as well. 
In this paper, we turn our attention also to Hahn and to Racah polynomials. The classification 
process becomes rather technical, however. Therefore, we have decided to present the proof of 
the complete classification only for dual Hahn polynomials (Section 3). For Hahn polynomials 
(Section 4) we give the final classification and corresponding two-diagonal matrices (but omit 
the proof), and for Racah polynomials we give the final classification and some examples of 
two-diagonal matrices in Appendix A. 

We should also note that the two-diagonal matrices appearing as a result of the doubling 
process are symmetric. So matrices like (1.3) do not appear directly but in their symmetrized 
form. Of course, as far as eigenvalues are concerned, this makes no difference (see Section 6 ). 

The doubling process of the polynomials considered here also gives rise to “new” sets of 
orthogonal polynomials. One could argue whether the term “new” is appropriate, since they 
arise by combining two known sets. The peculiar property is however that the combined set has 
a common unique weight function. Moreover, we shall see that the support set of these doubled 
polynomials is interesting, see the examples in Section 5. In this section, we also interpret 
the doubling process in the framework of Christoffel-Geronimus transforms. It will be clear 
that from our doubling process, one can deduce for which Christoffel parameter the Christoffel 
transform of a Hahn, dual Hahn or Racah polynomial is again a Hahn, dual Hahn or Racah 
polynomial with shifted parameters. 

In Section 6 we reconsider the two-diagonal matrices that have appeared in the previous 
sections. It should be clear that we get several classes of two-diagonal matrices (with parameters) 
for which the eigenvalues (and eigenvectors) have an explicit and rather simple form. This section 
reviews such matrices as new and potentially interesting examples of eigenvalue test matrices. 
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In Section 7 we explore relations with other structures. Recall that in finite-dimensional 
representations of the Lie algebra su(2), with common generators J+, J_ and Jo, the matrix 
of J + + J- also has a symmetric two-diagonal form. The new two-diagonal matrices appearing 
in this paper can be seen as representation matrices of deformations or extensions of su(2). We 
give the algebraic relations that follow from the “representation matrices” obtained here. The 
algebras are not studied in detail, but it is clear that they could be of interest on their own. 
The general algebras have two parameters, and we indicate how special cases with only one 
parameter are of importance for the construction of finite oscillator models. 


2 Introductory example 


We start our analysis by the explanation of a known example taken from [27]. For this example, 
we first recall the definition of Hahn and dual Hahn polynomials and some of the classical 
notations and properties. 

The Hahn polynomial Q n (x; a, j3, N ) [13, 19, 23] of degree n, n = 0,1,..., N, in the variable x, 
with parameters a > — 1 and (3 > — 1 (or a < —N and [3 < — N ) is defined by [13, 19, 23] 


Q n (x;a,P,N) 


3-^2 


/ —n, n + a + (3 + 1, — x 
V a + 1 ,—N 



( 2 . 1 ) 


Herein, the function 3 F 2 is the generalized hypergeometric series [4, 26] 


fa,b,c \ = (a)k{b)k(c) k z k 
V d,e ,Z ) ^ ( d) k (e) k k\' 


( 2 . 2 ) 


In (2.1), the series is terminating because of the appearance of the negative integer —n as 
a numerator parameter. Note that in (2.2) we use the common notation for Pochhammer 
symbols [4, 26] (a) k = a(a + 1) • • • (a + k — 1) for k = 1, 2,... and (a)o = 1. Hahn polynomials 
satisfy a (discrete) orthogonality relation [13, 19] 


N 


Y : w(x- a, /3, N)Q n (x ; a, /3, N)Q n >(x ; a, /3, N ) = h n (a, ft, N ) S n y, 


(2.3) 


x=0 


where 


w(x- a, 0, N) = ( Q ^ X ^j ( N + x X ), x = 0,1,..., N, 


h n (a,P,N) = 


(—l) n (n + a + /3 + 1)jv+i(/3 + 1 ) n n\ 


(2 n + a. + (3 + l)(a + l) n (—N) n N\ 
We denote the orthonormal Hahn functions as follows 


A , a , n \fw{x-, a, P, N) Q n (x ; a, /3, N ) 

Qn{ x j Oi, /3, N) = - — -. 

x/h n (a, (3, N) 

The Hahn polynomials satisfy the following recurrence relation [19, equation (9.5.3)] 
A (x)y n (x) = A(n)y n+ i(x) - (A(n) + C{n))y n (x) + C(n)y n -i(x) 


(2.4) 


Vn{x) = Q n (x;a,/3,N), A(x) = -x, 

. _ (n + a + l)(n + a + f3 + 1)(N - n) 
(2 n + 0 . + f3 + l)(2n + ot + [3 + 2) 


(2.5) 

n(n + a + (3 + N + l)(n + j3) 

(2 n + a + (3)(2n + a + (3 + 1) 


with 
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Related to the Hahn polynomials are the dual Hahn polynomials: R n (A(x); 7, 6, N) of deg¬ 
ree n, n = 0, 1,..., N, in the variable A(x) = x(x + 7 + 5 + 1), with parameters 7 > — 1 and 
5 > — 1 (or 7 < — N and 5 < — N ) which are defined similarly to (2.1) [13, 19, 23] 


R n (A(x); 7 ,5,lV) 


3-^2 


/— x, x + 7 + 5 + 1, —n 
v 7 + 1 ,-AT 



( 2 . 6 ) 


As is well known, the (discrete) orthogonality relation of the dual Hahn polynomials is just the 
“dual” of (2.3) 


N 


T: w(x; 7 , 5, N)R n (\(x)-, 7 ,5, N)R n ,(\(x)- 7 ,5, N) 

x=0 


^n('7 5 ^5 -^)^n,n / 5 


(2.7) 


where 


^0+7,5, AT) = 


/i n (7,5, AT) = 


(2,T + 7 + 5 + l )( 7 + l) :c (-Ar), T A^! 
(—l) 3 : (a: + 7 + 5 + 1)jv+i( 5 + l)+c! : 
7 + n\/CV + 5 — n 
n /V N — n 


1 -1 


Orthonormal dual Hahn functions are defined by 

5 r aa _ \M x ; 7 , < 5 , A^)R„(A(x); 7 , 5 , AT) 

R n {\(x)-'y, 6 ,N) = --■ 

K{ 7, <5, N) 


( 2 . 8 ) 


Dual Hahn polynomials also satisfy a recurrence relation of the form (2.4), with [19, equa¬ 
tion (9.6.3)] 

y n (x) = R n ( A(x); 7 , 5 , AT), A(x) = A(x) = x(x + 7 + 5 + 1), 

A(n) = (n + 7 + l)(n — N), C(n) = n(n — 5 — N — 1). (2.9) 


In [27], the following difference equations involving two sets of Hahn polynomials were derived 
(for convenience we use the notation Q n (x) = Q n (x ; ot,[ 3+1, N ) and Q n (x) = Q n (x; a+1, /3, AI)): 


(IV + p + 1 — x)Q n (x) - (IV - x)Q n (x + 1) = (n + Q: + 1)(/? 1 + /j + 1) Q r ,.(x), (2.10) 

a + 1 

(x + l)(9 n (x) - (a + x + 2)Qn(® + 1 ) = -(a + l)Q n (x + 1). (2.11) 


Writing out these difference equations for x = 0,1,..., A^, the resulting set of equations can 
easily be written in matrix form. For this matrix form, let us use the normalized version of the 
polynomials, and construct the following (2 N + 2) x (2N + 2) matrix U with elements 


U 2 x,N—n — U‘2x,N-\-n -\-1 — 


(- 1) 3 

V2 


■Q n (x\a,P + 1 ,N), 


U2x+1,N—n — U2x+l,N+n+l — 


:-i) a 

V2 


Qn(x; a + 1,13, N), 


( 2 . 12 ) 

(2.13) 


where x,n 6 {0,1,..., AT}. By construction, this matrix is orthogonal [27]: the fact that the 
columns of U are orthonormal follows from the orthogonality relation of the Hahn polynomials, 
and from the signs in the matrix U. Thus U T U = UU 1 = I, the identity matrix. 
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The normalized difference equations (2.10), (2.11) for x = 
matrix form. The coefficients in the left hand sides of (2.10), 
(2 N + 2) x (21V + 2)-matrix of the form 

/ 0 M 0 0 \ 

Mo 0 M\ 

M = 0 Mi 0 0 

M 2N 

y o m 2N o J 

with 


0,1,..., N can then be cast in 
( 2 . 11 ) give rise to a tridiagonal 


(2.14) 


M 2k = ^{k + a + l){N + f3 + l-k), M 2k+1 = y/(k + l)(N-k). (2.15) 


Suppose a > — 1, f3 > — 1 or a < —N — 1, j3 < —N — 1 and let U be the orthogonal matrix 
determined in (2.12), (2.13). Then [27] the columns of U are the eigenvectors of M, i.e., 

MU = UD , (2.16) 


where D is a diagonal matrix containing the eigenvalues of M 


D = diag(— cat, ..., -ei, -e 0 , e 0 , ei,..., e N ), 

e k = -\/(a + k + 1)(/3 + k + 1), k = 0,1,..., N. (2.17) 


Note that the eigenvalues of the matrix M are (up to a factor 2) the same as those of 
the matrix C 2 N+ 2 ( y a, /3), the two-parameter extension of the Sylvester-Kac matrix. As we will 
further discuss in Section 6 , the above result proves that the eigenvalues of C 2 N +2 (a, f3) are 
indeed given by (1.4). Even more: the orthonormal eigenvectors of M are just the columns 
of H. 

Another way of looking at (2.16) is in terms of the dual Hahn polynomials. Interchanging x 
and n in the expressions (2.12), (2.13), we have 


U 2 n t N—x — U 2n ^N+x+\ — 


(~l) f 

y/2 


R n (X(x)-a,/3 + 1, N), 


U 2 n+l,N—x — ^2n+l,iV4-x-|-l — 


HT 

V2 


R n (X(x)-a + 1 ,/3,N), 


(2.18) 

(2.19) 


where x,n £ {0,1,..., N}. In this way, each row of the matrix U consists of a dual Hahn 
polynomial of a certain degree, having different parameter values for even and odd rows. Now, 
the relation (2.16) can be interpreted as a three-term recurrence relation with M being the Jacobi 
matrix. Two sets of (dual) Hahn polynomials (with different parameters) are thus combined 
into a new set of polynomials such that the Jacobi matrix for this new set has a simple two- 
diagonal form, with simple eigenvalues. The pair of difference equations (2.10), (2.11) involving 
two sets of Hahn polynomials then corresponds to the following relations involving the dual 
Hahn polynomials R n (x ) = R n ( A(x); 7 , 5 + 1, N ) and R n {x ) = R n { A(.x); 7 + 1, 5, N ): 


(N + 5 + 1 - n)R n {x) -{N- n)R n+1 (x) = {x + 7 + ^. + 6 + ^ (s), (2.20) 

(7 + 1 ) 

(n + l)Rn(x) - (n + 7 + 2)R n+1 (x) = -(7 + l)R n+ i(x). (2.21) 


This is in fact a special case of the so-called Christoffel transform of a dual Hahn polynomial 
with its transformation parameter chosen specifically so that the result is again a dual Hahn 
polynomial (with different parameters). We will further elaborate on this in Section 5. 

This introductory example, taken from [27], opens the following question: in how many 
ways can two sets of (dual) Hahn polynomials be combined such that the Jacobi matrix is 
two-diagonal? This will be answered in the following section. 
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3 Doubling dual Hahn polynomials: classification 

The essential relation in the previous example is the existence of a pair of “recurrence rela¬ 
tions” (2.20), (2.21) intertwining two types of dual Hahn polynomials (or equivalently a couple 
of difference equations (2.10), (2.11) for two types of their duals, the Hahn polynomials). Let 
us therefore examine the existence of such relations in general. Say we have two types of dual 
Hahn polynomials with different parameter values for 7 and 5 (and possibly N ) denoted by 
R n (X(x)', r y,6,N) and R n (\(x); 7 , 5, N), that are related in the following manner 

a(n)R n (X(x)-'y,5, N) + b(n)R n+ i (A(x); 7, 6, N) = d(x)R n (X(x);^,5,N) (3.1) 

a{n)R n (X(x);j,5,N) + S(n)17 n+ i(A(x); 7, <5, N) = d(x)R n+ i(X(x);'y, 5, N). (3.2) 

If we want these relations to correspond to a matrix identity like (2.16), then it is indeed 
necessary that the (unknown) functions a(n), a(n), b(n) and b(n) are functions of n and not 
of x, and that d(x) and d(x) are functions of x and not of n. Of course, the parameters 7 , 5, N, 
7 , 5, N can appear in these functions. 

In order to lift this technique also to other polynomials than just the dual Hahn polynomials, 
say we have the following relations between two sets of orthogonal polynomials of the same class, 
denoted by y n and y n , but with different parameter values 

a{n)y n + b(n)y n+ i = d(x)y n , (3.3) 

a{n)y n + b(n)y n+ 1 = d(x)y n+1 , (3.4) 

where a, a, b, b are independent of x and d, d are independent of n. Although (3.3), (3.4) are not 
actual recurrence relations since they involve both y n and y n , we will refer to a couple of such 
relations intertwining two types of orthogonal polynomials as “a pair of recurrence relations”. 
When substituting (3.4) in (3.3), we arrive at the following recurrence relation for y n 

a{n) [a(n - l)y„_i + b(n - 1 )y n ] + b(n) [a(n)y n + b{n)y n+ 1 ] = d(x)d{x)y n . (3.5) 

In the same manner, y n can be eliminated to find a recurrence relation for y n 

a{n - 1) [a(n - l)y n -i + b(n - 1 )y n ] + b{n - 1) [a(n)y n + b(n)y n+ 1 ] = d(x)d{x)y n - (3.6) 

Of course, the orthogonal polynomials y n already satisfy a three-term recurrence relation of 
the form (2.4). A comparison of the coefficients of y n + 1 , y n -, Un-i in (3.5), (3.6) with the known 
coefficients given in (2.9) leads to the following set of requirements for a, a, 6 , 6 , d, d 


a(n)a(n — 1) = C(n), (3-7) 

a(n — 1 )a(n — 1) = C(n), (3.8) 

a(n)b(n — 1) + a(n)b(n) — d(x) d{x) = — [A(x) + A(n) + C(n)] , (3.9) 

a{n)b{n — 1) + a(n — 1 )b(n — 1 ) — d{x)d(x) = — [A(x) + A(n) + C(ri)] , (3.10) 

b(n)b(n) = A(n), (3.11) 

b(n)b(n — 1) = A(n). (3-12) 

After a slight rearrangement of terms in the requirements (3.9) and (3.10), we arrive at two new 
equations where the left hand side is independent of x while the right hand side is independent 
of n, namely, 

a(n)b(n — 1 ) + a(n)b(n) + A(n) + C{n) = d{x)d{x) — A(x), (3.13) 

a(n)b(n — 1 ) + a(n — 1 )b(n — 1) + A{n) + C(n) = d(x)d{x) — A(x). (3-14) 
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Hence, the two sides must be independent of both n and x. By means of (3.7)-(3.12) we can 
eliminate A. A, C , C to find 

a(n) [a(n — 1) + b(n — 1)] + b(n) [a(n) + 6(n)] = d(x)d(x) — A(x), 

a(n — 1) [a(n — 1) + b(n — 1)] + b(n — 1) [a(n) + b(n)] = d{x)d(x) — A(x). 

Moreover, subtracting one from the other yields 

A(x) — A(x’) = a(n — 1) [a(n) — a(n — 1) — b(n — 1)] + b{n) [a(n) + b(n ) — b(n — 1)]. (3.15) 

Now, for a given class of orthogonal polynomials with recurrence relation of the form (2.4), 
we determine all possible functions a, a, b, b, d, d satisfying the list of requirements (3.7)-(3.12). 
Hereto, we proceed as follows 

• From (3.7) and (3.8) we observe that, up to a multiplicative factor, C(n ) is split into two 
functions, a(n — 1) and a(n — 1). When a(n — 1) is shifted by 1 in n and multiplied again 
by a(n — 1) we must arrive at C(n). Hence, C and C consist of an identical part, and 
a part which differs by a shift of 1 in n. This observation gives a first list of possibilities 
for a and a. 

• Similarly we find a list for b and b by means of (3.11) and (3.12). 

• These possibilities are then to be compared with requirements (3.9) and (3.10). From 
(3.13), (3.14) and (3.15) we get an expression for the product d(x)d(x). Finally, the set 
of remaining choices for a, a, 6, b are to be plugged in (3.4) and (3.3) in order to get d, d 
and to verify if these relations indeed hold. 

The actual performance of the procedure just described is still quite long and tedious, when 
carried out for a fixed class of polynomials. In what follows we achieve this for the dual Hahn 
polynomials, which have the easiest recurrence relation, and it takes about three pages to present 
this. The reader who wishes to skip the details can advance to Theorem 1. 

For dual Hahn polynomials, the data is given by (2.9) 

Un = Rn(A(x); 7 , 6, N), y n = R n { A(x); 7, <5, IV), A(x) = X(x) = x(x + 7 + 5 + 1), 
A(n ) = (n + 7 + l)(n — N), C (n) = n{n — 6 — N — 1), 

and with similar expressions for A(x), A(n) and C(n ) (with x, 7 , 6, N replaced by x, 7 , 5, N ). 
From (3.15), the following expression must be independent of x 

A(x) — A(x) = x(x + 7 + d + 1) — x(x + 7 + <5 + l). 

In order for the term in x 2 to disappear, we must have x = x + £ which gives 

x(x + 7 + 5 + 1 ) — (x + £)(x + £ + 7 + 5 + 1) = (7 + <5 — 7 — <5 — 2 £)x — £(£ + 7 + <5 + 1) 

and as we require the coefficient of x to be zero we find the following condition for £ 


7 + d- (j + S) = 2£. (3.16) 

From (3.8) we see that we have four distinct possible combinations for a(n — 1) and a(n — 1) 

a(n — 1) = lc a , a(n — 1) = n(n — 5 — N — l)c“ 1 , (al) 

a(n — 1) = nc a , a(n — 1) = (n — <5 — N — l)c“ 1 , (a2) 

a(n — 1 ) = (n — 5 — N — l)c a , a(n — 1) = nc~ l , (a3) 

a(n — 1) = n(n — 6 — N — 1 )c a , a(n — 1) = lc” 1 , (a4) 
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with c a a factor. Combining this with (3.7) we must have 
a(n)a{n — 1) = C(n) = n(n — 8 — N — l). 


This immediately implies that c a is independent of n, and (al)-(a4) yield the following possi¬ 
bilities 


n(n — 5 — N — 1) = n(n — 5 — N — 1) 

(n + l)(n — 5 — N — 1) = n(n — 5 — N — 1) 
(n — 5 — N)n = n(n — 5 — N — 1) 

(n + l)(n — 5 — N) = n(n — 6 — N — 1) 


5 + N = 5 + N, (aT) 

S + N + 1 = 0 A 5 + N- 1-2 = 0, (a2 7 ) 

5 + N = S + N + 1, (a3') 

6 + N = 0 A 6 + N + 2 = 0. (a4 7 ) 


Because of the restriction on 5 the option (a4 7 ) is ineligible, leaving (al / )-(a3 / ) as only viable 
options. 

In a similar way, from (3.12) we see that we have four possible combinations for b(n ) and b(n), 


b(n) 

= 1 c b , 

b(n - 1) = (n + 7 + l)(n - N)c b 

(bl) 

b(n ) 

= (n + 7 + l)c 6 , 

b(n — 1) = (n — N)c b l , 

(b 2 ) 

b(n) 

= (n - N)c b , 

b(n- 1) = (n + 7 + l)c^ 1 , 

(b3) 

b(n) 

= (n + 7 + l)(n — N)c b , 

b(n~ 1) = lc b 1 - 

(b4) 


Combining this with (3.11) we must have 

b(n)b(n) = A(n) = (n + 7 + l)(n — N) 

This implies that q, is independent of n and moreover for (bl)-(b4) yields 

(n + 7 + 2 )(n — N + 1) = (n + 7 + l)(n — N) 

(n + 7 + l)(n — N + 1) = (n + 7 + l)(n — N) 

(n + 7 + 2 )(n — IV) = (n + 7 + l)(n — IV) 

(n + 7 + l)(n — N) = (n + 7 + l)(n — N) 

We thus have four viable options for b, b and three for a, a, giving a total of 12 possible 
combinations, which we will systematically consider and treat. 

Case (bl). Plugging (bl) in (3.14), we get 

a(n)(n + 7 + l)(n — N) c^ 1 + a(n — l)c& + (n + 7 + l)(n — N) + n(n — 5 — N — 1) 

= d(x)d(x) — A(x). 


> 7 + 1 = 7 A JV- 1 = iV, (bl') 

> 7 = 7 A A - 1 = IV, (b2') 

> j + 1= jAJV = JV, (b3') 

> 7 = 7 a N = N. (b4') 


As the right hand side is independent of n, so must be the left hand side. This eliminates op¬ 
tions (a2) and (a3) for a, a as that would result in a third order term in n which cannot vanish. 
On the other hand, (al) yields 

(n + 7 + 1 )(n — A^)— + n(n — <5 — N — 1)— + (n + 7 + l)(n — N) + n(n — 5 — N — 1) 

Cb Ca 

= d{x)d{x) — A(x). 

This must be independent of n, so the coefficient of n 2 in the left hand side must vanish, hence 
c a /cb + c-b/ca + 2 = 0 or thus c a /cfe = —1. For this value of c a /c& the left hand side equals 
zero and is indeed independent of n. Note that this leaves one degree of freedom as only the 
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ratio c a /cb is fixed. This is just a global scalar factor for (3.3) and (3.4), also present in (2.4). 
Henceforth, for convenience, we set c a = 1 and q, = —1. 

The combined options (bl) and (al) thus give a valid set of equations of the form (3.3) 
and (3.4), and they correspond to the parameter values 

7 = 7 + 1, <5 = 6 + 1, N = N — 1. 

Moreover, by means of (3.16) we find £ = — 1 and so x = x — 1. Finally, plugging these a, a, 
b, b in (3.3) and (3.4), and putting n = 0 we find 

i?o(A(x); 7 , <5, N) - R 1 (A(s); 7 ,5, N) = ^^ = d(x) 

and similarly d(x) = V(y + 1). Hence, for R n (x) = R n { A(x); 7 , 6, N ) and R n (x) = R n (X(x — 1); 
7 + 1 , 5 + 1, N — 1) we have the relations 

R n \x) R n +i[x) — N(^y + 1) 

-(n + 1)(N - n + 6)R n {x) + (N - n - 1 )(n + 7 + 2 )R n+i (x) = ^(7 + l)R n+ i(x). 

Interchanging x and n, these recurrence relations for dual Hahn polynomials are precisely the 
known actions of the forward and backward shift operator for Hahn polynomials [19, equa¬ 
tions (9.5.6) and (9.5.8)]. 

Case (b2). Next, we consider the option (b2) for b , b. Plugging (b2) in (3.14), we get 

a(n)(n — N)^ 1 + a{n — l)(n + 7 ) 0 , + (n + 7 + l)(n — N) + n(n — 5 — N — 1) 

= d{x)d{x) — A(x). 

Since the left hand side must be independent of n, option (al) is ruled out. Also option (a 2 ) 
is ruled out: using (a2) and <5 + N + 1 = 0 (from (a2')), the left hand side again cannot be 
independent of n. Only (a3) remains, giving 

(n — 5 — N)(n — N) — + n(n + 7 )— + (n + 7 + l)(n — N) + n{n — 5 — N — 1) 

Cb Ca 

= d(x)d(x) — A(x). 

In order for n 2 in the left hand side to vanish, we again require c a /cb = —1. This gives 
—N(N + 7 + d + 1) = d{x)d{x) — A(x), 


and we see that both sides are indeed independent of n. 

The combined options (b2) and (a3) also give a valid set of equations of the form (3.3) 
and (3.4), now corresponding to the parameter values 

^ = 7j 6 = 6, N = N- 1. 


Moreover, by means of (3.16) we find ^ = 0 and so x = x. Putting again n = 0 in (3.3) and (3.4) 
for these a, a, b, b we find 

(-6 - N)R 0 (A(x);j, 6, N) - (7 + l)Ri(X(x)\ 7, 6, N ) 

(N - x)(x + 7 + 6 + N + 1 ) 3 

= -= d{x) 


N 
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and similarly d(x) = N. The relations in question are then, for R n (x) = i? n (A(x); 7 ,5, N) and 
R n {x) = Rn(X(x)-,j,6,N - 1) 

(' n-8 - N)R n (x) - (n + 7 + l)R n+1 (x) = + S + N + 

(n + 1 )R n {x) - (n - N + l)R n+ i(x) = NR n+1 (x), 

which can be verified algebraically or by means of a computer algebra package. 

Case (b3). The next option to consider is (b3), for which (3.14) becomes 

a(n){n + 7 + l )^ 1 + a(n — l)(n — N — l)q, + (n + 7 + l)(n — N) + n{n — 5 — N — 1) 

= d{x)d{x) — A(x). 

The independence of n in the left hand side again rules out options (al) and (a2), while (a3) 
gives 


{n — 5 — N){n + 7 + 1)— + n{n — N — 1)— + {n + 7 + l)(n — N ) + n{n — S — N — 1) 

Cfr c a 

= d(x)d(x) — A(x). 

Also here, we require c a /c\, = — 1 to arrive at a left hand side independent of n, namely 
(7 + 1)5 = d(x)d(x) — A(x). 

The combined options (b3) and (a3) thus give a valid set of equations of the form (3.3) and (3.4), 
and they correspond to the parameter values 

7 = 7 + 1 , <5 = 5 — 1 , N = N] 

by means of (3.16) we find ^ = 0 and so x = x. Finally, plugging these a, a, 6 , b in (3.3) 
and (3.4) and putting n = 0 we find 

(S ~ N)R 0 (X(xy, 7 , 5, N ) + NRi(X(x);j, 5, N ) = _ (a; + 7 + 1 )(* + S) = ^ 

(7 + 1 ) 

and similarly d(x) = 7 + 1 . 

Hence we have the relations, for R n (x) = R n (X(x); r y,6,N) and R n {x) = R n (X(x);'y + 1, 

5 -1,N) 

-( n-S- N)R n (x ) + (n - N)R n+1 (x) = ^ + 7 . + ^ + ^ R n (x), 

(7 + 1) 

~(n + 1 )Rn(x) + (n + 7 + 2)R n+ i(x) = (7 + l)R n+1 (x). 

These can again be verified algebraically or by means of a computer algebra package. Note that 
these relations coincide with (2.20), (2.21) from the previous section (up to a shift 5 —> 5 + 1). 

Case (b4). The final option (b4) for b , b does not correspond to a valid set of equations of 
the form (3.3) and (3.4) as the left hand side of (3.14) can never be independent of n for either 
options (al), (a2) or (a3). 

This completes the analysis in the case of dual Hahn polynomials, and we have the following 
result 
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Theorem 1. The only way to double dual Hahn polynomials, i.e., to combine two sets of dual 
Hahn polynomials such that they satisfy a pair of recurrence relations of the form (3.1), (3.2) is 
one of the three cases: 

dual Hahn I, R n (x) = R n ( A(x); 7 , 5, N ) and R n {x ) = R n (\(x — 1) 57 + 1 ,5 + 1,IV — 1): 

Rn(X) Rn+1\X) — N[pf + 1 ) ^n\X), 

-{n + 1)(N - n + 5)R n (x ) + (N - n - 1 )(n + 7 + 2)R n+ i{x) = N( 7 + l)R n+ i(x). 
dual Hahn II. R n {x) = _R n (A(x); 7, 5, IV) and R n (x) = 7? n (A(x); 7, 5, IV — 1); 

(n — 5 — N)R n (x) - (n + 7 + l)i? n+ i(x) = - —— + ^ ^ 

(n + l)R n (x) — (n - IV + l)R n+ i(x) = lVR n+ i(x). 
dual Hahn III, R n (x) = R n { A(x); 7 , 5, N ) and R n (x) = R n (X(x)] 7 + 1,5—1, IV): 

-(n - 5 - N)R n (x) + (n - lV)R n+ 1 (x) = (X + 7 . + ^ + ^ R n (^), 

(7 + 1) 

-(n + l)R„(x) + (n + 7 + 2 )R n+ i(x) = (7 + l)R n+ i(x). 


By interchanging x and n, each of the recurrence relations for dual Hahn polynomials in the 
previous theorem gives rise to a set of forward and backward shift operators for regular Hahn 
polynomials. The case dual Hahn I corresponds just to the known forward and backward shift 
operators for Hahn polynomials [19]: Q n (x ) = Q n (x; a, /3, N) and Q n (x ) = Q n (x\a + l,/3 + 1, 
IV- 1): 


Qn{x ) Qn(x + 1 ) — 


n(n + a + j3 + 1) 
JV(a + l) 


Qn— 1 (®)) 


-(x + 1)(1V - X + /3)Qn-i(®) + (iv — X — l)(x + a + 2)(5 rt -i(£ + 1) 
= lV(a + 1 )Qn(x + 1). 


The case dual Hahn III corresponds to our introductory example (2.10), (2.11) (up to a shift 
(5 —> j3 + 1), and appears already in [27]. The case dual Hahn II yields a new set of relations 
(encountered recently in [16, equations (16), (17)]), namely Q n {x ) = Q n (x; a , /3, N ) and Q n {x ) = 
Q n (x;a,/3,N - 1 ): 

(x - P - N)Q n (x ) - (x + a + l)Q n (a; + 1) = — n )( n + ( * + P + N + 1 - > Q n (x), 

(x + 1 )Qn(x) - (x - N + 1 )Qn(x + 1) = NQ n (x + 1). 


The most important thing is, however, that we have classified the possible cases. 

Because the sets of recurrence relations are of the form (3.1), (3.2), they can be cast in matrix 
form, like in (2.16), with a simple two-diagonal matrix. For the case dual Hahn I, note that 
the IV-values of R n (x) and R n (x ) differ by 1, so the definition of the matrix U (again in terms 
of the normalized version of the polynomials) requires a little bit more attention. The matrix U 
is now of order ( 2 N + 1 ) x ( 2 N + 1 ) with matrix elements 


U2n,N—x — U2n,N+x — 


(~l) r 

v /2 


R n ( A(x); 7 ,5,1V), 


x = 1,... ,77, 




2n+l,N—x 


-x = -u, 


2n-\-l,N-\-x 


(- 1 )’ 

V 2 


-R n ( A(x - 1); 7 + 1,5 + 1,1V - 1), x = 1,... ,1V, 


U 2 n,N = (—l) n -R n (A(0); 7 ,5, IV), U 2 U +i,n = 0, 


(3.17) 
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where the row index of the matrix U (denoted here by 2n or 2n + 1, depending on the parity 
of the index) also runs over the integers from 0 up to 2N. This matrix U is orthogonal: the 
orthogonality relation of the dual Hahn polynomials (2.7) and the signs in the matrix U imply 
that its rows are orthonormal. Thus U T U = UU T = /, the identity matrix. Then the recurrence 
relations for dual Hahn I of Theorem 1 are now reformulated in terms of a two-diagonal 
(2 N + 1) x (21V + l)-matrix of the form 


/ 0 M 0 0 \ 

Mo 0 Mi 

M = 0 Mi 0 0 

M 2 7V-1 

\ 0 M 2 n- i 0 ) 

Explicitly 


(3.18) 


Proposition 2 (dual Hahn I). Suppose 7 > —1, 5 > —1. Let M be the two-diagonal mat¬ 
rix (3.18) with 

M 2k = v / (fc + 7 + l)(Ar_fe), M 2k+ 1 = VC k + l){N + S-k ), (3.19) 


and U the orthogonal matrix determined in (3.17). Then the columns of U are the eigenvectors 
of M, i.e., MU = UD, where D is a diagonal matrix containing the eigenvalues of M 


D = diag(— cat, ..., -ei,0,ei,.. -,e N ), 

e k = k(k -\- 'y 5 -\- 1), k = 1,..., N. (3.20) 


Note that we have kept only the conditions under which the matrix M is real. The other 
conditions for which the dual Hahn polynomials in (3.17) can be normalized (namely 7 < —N, 
5 < —N ) would give rise to imaginary values in (3.19). In such a case, the relation MU = UD 
remains valid, and also D would have imaginary values. 

For the case dual Hahn II, the matrix U is again of order (2 N + 1) x (2 N + 1) with matrix 
elements 


U 2n , x = U 2n , 2N - x = -j=Rn(\(x); 7 , 5,N), x = 0,..., N — 1, 


U 2 n+l,x — U2n+l,2N—x — 1 ), 

U 2n> N = Rn(^(N); 7 , S, N), U 2 n -\.x,N = 0, 


x = 0,..., N — 1, 


(3.21) 


where the row indices are as in (3.17). The orthogonality relation of the dual Hahn polynomials 
and the signs in the matrix U imply that its rows are orthonormal, so U T U = UU T = I. The 
pair of recurrence relations for dual Hahn II of Theorem 1 yield 

Proposition 3 (dual Hahn II). Suppose 7 > —l, 5 > — 1. Let M be a tridiagonal (2 N + 1) x 
(2 N + l)-matrix of the form (3.18) with 

M 2k = y/(N + 6-k)(N-k ), M 2k+l = v / (fc + l)(fe + 7 + l), (3.22) 


and U the orthogonal matrix determined in (3.21). Then the columns of U are the eigenvectors 
of M, i.e., MU = UD, where D is a diagonal matrix containing the eigenvalues of M 


D = diag(— cat, ..., -ei, 0 , ei,..., cat), 
e k = yjk{pi + 5 + 1 + 2 N — k ), 


k = 1,... ,N. 
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Note that the order in which the normalized dual Hahn polynomials appear in the matrix U 
is different for (3.17) and (3.21). This is related to the indices of the polynomials in the relations 
of Theorem 1. 

Finally, for the case dual Hahn III, the matrix U is given by (2.18), (2.19) and we reca¬ 
pitulate the results given at the end of the previous section, now in terms of the dual Hahn 
parameters 7 and 5. 

Proposition 4 (dual Hahn III). Suppose 7 > —1, h > —1 or 7 < —N — 1, 5 < —N — 1. Let M 
be the tridiagonal matrix (2.14) with 

M 2k = ^{k + 1 + l){N + 5 + l-k), M 2k+ 1 = sj{h + l)(N-k), (3.23) 

and U the orthogonal matrix determined in (2.18), (2.19). Then the columns of U are the 
eigenvectors of M, i.e., MU = UD, where D is a diagonal matrix containing the eigenvalues 
of M 


D = diag(—ejv, - - -, —ei, -eo, e 0 , ei, •.., e N ), 

= 's/ (k + 7 + l)(fc + <5 + 1), k = 0,l,...,N. 

To conclude for dual Hahn polynomials: there are three sets of recurrence relations of the 
form (3.1), (3.2). Each of the three cases gives rise to a two-diagonal matrix with simple and 
explicit eigenvalues, and eigenvectors given in terms of two sets of dual Hahn polynomials. 


4 Doubling Hahn polynomials 


The technique presented in the previous section can be applied to other types of discrete orthogo¬ 
nal polynomials with a finite spectrum. We have done this for Hahn polynomials. One level up 
in the hierarchy of orthogonal polynomials of hypergeometric type are the Racah polynomials. 
Also for Racah polynomials we have applied the technique, but here the description of the results 
becomes very technical. So we shall leave the results for Racah polynomials for Appendix A. 

For Hahn polynomials the analysis is again straightforward but tedious, so let us skip the 
details of the computation and present just the final outcome here. Applying the technique 
described in (3.3)—(3.15), with y n = Q n (x\ a, (3, N) and y n = Q n {x;a, f3, N) yields the following 
result. 


Theorem 5. The only way to combine two sets of Hahn polynomials such that they satisfy a pair 
of recurrence relations of the form (3.3), (3.4) is one of the four cases: 

Hahn I, Q n {x) = Q n (x;a, /3, N) and Q n {x ) = Q n (x;a + l,/3,N): 


(n + a + (3 + N + 2) 
{fin + ol + /3 + 2 ) 
(n + l)(n + /3 + 1) 
(2 n + cy. + /3 + 3) 


~Qn(x ) 
Qn(x ) + 


(N-n) „ 

( 2 n + a + /5 + 2 )^ n+lW_ 
(n + a + /3 + 2)(n + a + 2) 
(2n + ol + f3 + 3) 


(a + x + 1 ) 
(a + 1 ) 


Qn {x) 


Q n +i{x) = (a + l)Q n+ i(x). 


Hahn II, Q n (x) = Q n (x;a,/3,N) and Q n {x) = Q n (x - 1; a + 1, /3, N - 1): 

1 „ , . 1 , , x 


Qn{x) 


{In + 01 + /3 + 2) {In + ex + /3 + 2) 

(n + l)(n + P + l)(n + a + + N + 2) 


Qn+l(x ) = 


N{a + 1) 


Qn{x ), 


+ 


(2n + a + (3 + 3) 

(n + a + /3 + 2)(N — n — l)(n + a + 2) 
(2 n + a + (3 + 3) 


Qn+i{x) = N{a + l)Q n+ i(x). 
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Hahn III, Q n (x) = Q n {x\ a, f3, N) and Q n (x) = Q n (x-,a,/3 + 1,1V): 

(n + P + l)(n + N + 2 + at + P) ^ f ^ | (IV — n)(n -\- ot -\- Y) ^ ^ ^ 

^cn\X) + , 0 | : 7> : 777 ^Jn+l\X ) 


(2 n + a + P 4- 2) 


(2 n + ol + P + 2) 

= (P + 1 + N - x)Qn(x), 

(n + 1) A/'\ (ra + cc + /3 + 2) ~ 

-Qn(®) + 77:-;- , 0 , L Qn + l(x) = Qn+l(x) 


(2 n + a + P + 3) (2 n + a + P + 3) 

Hahn IV, Q n (x) = Q n (x-,a,P,N ) and Q n (x) = Q n (x;a,P + 1, N - 1): 

(n + ' 9 + 1) -o„(x) + <" +a tV .» <?»+■(*) - HH.W. 


( 2 n + ol + P + 2 ) 


(2 n + ol + P + 2 ) 


V 


(?r + l)(n + a +/3 + IV + 2) ^ ^ , (IV — n — l)(n + a + P + 2)^ t„\ _ AT n i„\ 

^Jn(x) -p . , Qn+l(2') — dy Cjn-\-l\X) 


(2n + ol 4- P + 3) 


(2 n + a + P + 3) 


Note that when interchanging x and n the relations in Hahn II coincide with the known 
forward and backward shift operator relations for dual Hahn polynomials [19, equations (9.6.6) 
and (9.6.8)]. In the same way, the other cases yield new forward and backward shift operator 
relations for dual Hahn polynomials. 

Since the recurrence relations are of the form (3.3), (3.4), they can be cast in matrix form 
with a two-diagonal matrix. We shall write the matrix elements again in terms of normalized 
polynomials. For the case Hahn I, the matrix U of order (2 N + 2) x (21V + 2), with elements 


U2n,N—x — U r 2n.N-\-x -\-1 — 


(~l) r 

V 2 


Q n (x;a,P,N), 


U2n+l,N—x — ^2re+l,iV+®+l — 


:-D r 

V2 


Qn(x', OL + 1 ,P,N) 


(4.1) 


where x, n € {0,1,..., IV}, is orthogonal, and the recurrence relations yield 


Proposition 6 (Hahn I). Suppose that 7 ,6> —1. Let M be a tridiagonal (2 N + 2) x (2 N + 2)- 
matrix of the form (2.14) with 


M'2k — 


(k + a + l)(k + a + P + l)(k + a + P + 2 + N ) 
(2 k + a + P + l)(2fc + q + P + 2) 


M'2k+\ — 


1 (k + p + l)(k + l)(N-k) 

(2k + ol + p + 2)(2 k + 0. + P + 3) 


and U the orthogonal matrix determined in (4.1). Then the columns of U are the eigenvectors 
of M, i.e., MU = UD, where D is a diagonal matrix containing the eigenvalues of M 

D = diag(—ejv,..., -ei, -e 0 , £ 0 , 6 !,..., e N ), 

ek = Vk + a + 1 , k = 0 , 1 , ..., N. (4.2) 

For the case Hahn II, the orthogonal matrix U is of order (2 N +1) x (2 N + 1), with elements 

(-1 r 


U 2 n,N—x — U2n,N+x — 


V2, 


Qn(x] a, P, N), 


x = 1,... ,1V, 


U2n+l,N—x — Lf2n+l,N-\-x — 


nr 

V 2 


Qn(x - l;a + 1,P, N - 1), x = l , ...,1V, (4.3) 


U2n,N = (~l) n Qn(0;a,P,N), U2n+1,N = 0, 

where the row indices are as in (3.17). The recurrence relations for Hahn II yield 
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Proposition 7 (Hahn II). Suppose that a,/3 > —1 or a, (3 < —N. Let M be a tridiagonal 
(2 N + 1) X (2N + 1)-matrix of the form (3.18) with 



(k + a + 1) (A: + a + f3 + 1)(IV — k ) 

(2 k + a + /3 + 1)(2 k + a + /3 + 2) ' 


(fc + /3 + l)(fe + o; + /3 + 2 + Ll)(k + 1) 
(2k + a + f3 + 2)(2k + a + /3 + 3) : 


and U the orthogonal matrix determined in (4.3). Then the columns of U are the eigenvectors 
of M, i.e., MU = UD, where D is a diagonal matrix containing the eigenvalues of M: 

.D = diag(-ejv,. • • ,-ei,0,ei,... ,ejv), e k = Vk , k = l,...,N. 


Note that for both cases, the two-diagonal matrix M becomes more complicated compared 
to the cases for dual Hahn polynomials, but the matrix D of eigenvalues becomes simpler. 

For the two remaining cases we need not give all details: the matrix M for the case Hahn III 
is equal to the matrix M for the case Hahn I with the replacement af>j3, and so its eigenvalues 
are ± \/k + /3 + 1, k = 0,1,..., N. And the matrix M for the case Hahn IV is equal to the 
matrix M for the case Hahn II with the same replacement a •<->•/?, so its eigenvalues are 0 
and P'/k, k = 1,..., N. 


5 Polynomial systems, Christoffel and Geronimus transforms 

So far, we have only partially explained why the technique in the previous sections is referred to as 
“doubling” polynomials. It is indeed a fact that the combination of two sets of polynomials, each 
with different parameters, yields a new set of orthogonal polynomials. This can be compared to 
the well known situation of combining two sets of generalized Laguerre polynomials (both with 
different parameters a and a — 1) into the set of “generalized Hermite polynomials” [8]. There, 
for a > 0, one defines 


P% 


' (X) = \lWn L ^ 1>{x2) ' P2U+l{X) = \iuU 1 XL( “ >{x2) - 


(5.1) 


Then the orthogonality relation of Laguerre polynomials leads to the orthogonality of the poly¬ 
nomials (5.1): 


r+oO 


w(x)P n (x)P n '(x)dx = T(a)d„ 


where 


w(x) = e- x2 \x\ 2a ~ 1 . (5.2) 

Note that the even polynomials are Laguerre polynomials in x 2 (for parameter a — 1), and the 
odd polynomials are Laguerre polynomials in x 2 (for parameter a) multiplied by a factor x. 
The weight function (5.2) is common for both types of polynomials. It is this phenomenon that 
appears here too in our doubling process of Hahn or dual Hahn polynomials. 

From a more general point of view, this fits in the context of obtaining a new family of 
orthogonal polynomials starting from a set of orthogonal polynomials and its kernel partner 
related by a Christoffel transform [8, 22, 32]. In a way, our classification determines for which 
Christoffel parameter v (see [32] for the notation) the Christoffel transform of a Hahn, dual 
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Hahn or Racah polynomial is again a Hahn, dual Hahn or Racah polynomial with possibly 
different parameters. This determines moreover quite explicitly the common weight function. 

For a dual Hahn polynomial R n (x) = R n (X(x); 7 , 5, N), with data given in (2.9), and a Chris- 
toffel parameter v the kernel partner is given by the transform 

D (^\ Rn+i(x) — a n R n (x ) ^ R n+ \{v) 

n\ x ) — 7"7 \ 77 7 ; a n ~ D / \ • l^.oj 

A(x) - A(u) R n {v) 

Because of the recurrence relation (2.4) and what is called the Geronimus transform the original 
polynomials can also be expressed in terms of the kernel partners. This is usually done for 
monic polynomials (see [32, equations (3.2) and (3.3)]), but it can be extended to non-monic 
dual Hahn polynomials as follows 

R n (x) = A(n)P n (x) - b n P n -i(x) (5.4) 

where the coefficients b n are related to the recurrence relation (2.4) as follows 

b n a n -i = C(n), A(n)a n + b n = A{n) + C(n) + A(u). (5.5) 


Our classification now shows that only for v equal to one of the values 0, N or — 6 , the kernel 
partner P n {x) will again be a dual Hahn polynomial. Indeed, taking for example v = 0 in (5.3) 
we have R n ( 0) = 1 and 


P n {x) 


Rn+\(x) Rn(x') 

A(x) 


i^(A(x - 1); 7 + M + 1, N- 1 ), 


where we used the first relation of dual Hahn I to obtain again a dual Hahn polynomial. The 
reverse transform (5.4) follows immediately from the second relation of dual Hahn I. Similarly, 
taking v = N in (5.3) we have R n (N) = (-N — S) n / (7 + l) n and 


R n+ i(x) - Rn(x)(n - 6 - N)/(n + 7 + 1) 
n[X ~ (x- N)(x + N + 7 + 6 + 1) 


-1 

N(n + 7 + I) 


Rn{ X{x);n/,S,N - 1 ), 


which we obtained using the first relation of dual Hahn II. For the reverse transform (5.4) we 
find, using the second relation of dual Hahn II with shifted n i— > n — 1, 


A(n)P n (x) - b n P n - i(x) 


~jy fi n(A(x); 7,(5, - 1) + -^R n _i(A(x); 7 , <5, N - 1) 

R n {x). 


For the last case, taking v = —5 in (5.3) we have R n (—5) = (-N — 5) n /(—N) n and 
R n +i(x) - R n (x)(n - 5 - N)/(n - N) 


P n (x) = 


1 


(x + 7 + l)(x + (5) 

-Rn( A(x );7 + 1,(5- 1, N), 


(7 + l)(n - N )' 

which we obtained using the first relation of dual Hahn III. For the transform (5.4) we have 

ti 'y “l - 1 

A(n)P n (x) - b n P n -i(x) = -—— Rn( A(x );7 + 1,(5 - 1, N - 1) 


7+1 

n 

7 + 1 


R„_i(A(x );7 + 1,(5- 1,1V), 


which equals R n {x) by the second relation of dual Hahn III. 
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In a similar way, for the Hahn polynomials, putting Q n (x) = Q n (x; a,b, N), using the 
data (2.5) in 

j-) f ^ Qn+l(%) 0 ’nQn(%) _ Qn+ 1 (^) 

n[X) ~ A(s)-AW ’ Qn{v) ’ 

and in (5.4), (5.5), the cases Hahn I, II, III, V correspond respectively to the choices —a — 1, 
0, IV + P + 1 and N for v. 

The task of determining for which Christoffel parameter v the kernel partner of a dual Hahn 
polynomial is again of the same family is not trivial. It comes down to finding a pair of recurrence 
relations of the form (3.1), (3.2) with coefficients related to v as in (5.3). We have classified 
these for general coefficients, without a relation to v, and we observe that each solution indeed 
corresponds to a specific choice for v. 

The transforms (5.3), (5.4) give rise to new orthogonal systems, but in general there is no way 
of writing the common weight function. However, since here both sets are of the same family, 
we can actually do this. Let us begin with the dual Hahn polynomials, in particular the case 
dual Hahn I, for which the corresponding matrix U is given in (3.17). They give rise to a new 
family of discrete orthogonal polynomials with the relation MU = UD corresponding to their 
three term recurrence relation with Jacobi matrix M (3.19). In general the support of the weight 
function is equal to the spectrum of the Jacobi matrix [5, 18, 20, 21], After simplifying with the 
normalization factors ( 2 . 8 ), this leads to a discrete orthogonality of polynomials, with support 
equal to the eigenvalues of M (so in this case, the support follows from (3.20)). Concretely, for 
the case under consideration, we have 


Proposition 8 . Let 7 > —1, 5 > —1, and consider the 2 N + 1 polynomials 
P 2 n{q) = —n = 0, 1,. ..,N, 


P 


2n+l 


v/2 


n = 0 , 1 , 


,1V- 1. 


These polynomials satisfy the discrete orthogonality relation 


E 


(-l ) fc (2 k + 7 + & + 1)(7 + 1 )k(~N) k N\ 
(k + 7 + 6 + 1)at+i(<5 + 1 )kkl 


(1 + 5 q ,o)P n (q)P n '(q) 


h+[n/2\\f5 + N-[n/2\ 
V Ln/2j J\ N — \n/2\ 



5 


n,n' 


(5.6) 


with 


S= {0,±^k(k+ 1 + 5 + 1), k = 1 , 2 ,. ..,1V}. 

Note that for q E S, q 2 = k(k + 7 + <5 + 1) = A (k), and the polynomial P 2 n{q) is of the form 
R n {\(k)\ 7, 5, N ). In that case, q 2 — 7 — 6 — 2 = (k — l)({k — 1) + (7 +1) + (<5 +1) +1) = A (k — 1), 
so the polynomial P 2 n+i(^) is of the form R n (X(k — 1);7+1,J + 1,1V — 1). The interpretation 
of the weight function in the left hand side of (5.6) is as follows: each q in the support S is 
mapped to a lc-value belonging to {0,1,..., IV}, and then the weight depends on this k- value. 

Now we turn to the classification of Section 4, where the corresponding orthogonal matri¬ 
ces U are given in terms of (normalized) Hahn polynomials, for the case Hahn I, the matrix U 
is given in (4.1), and the spectrum of the matrix M is given by (4.2). After simplifying the 
normalization factors, the orthogonality of the rows of U gives rise to 
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Proposition 9. Let a > — 1, f3 > — 1, and consider the 2N + 2 polynomials, n = 0,1,..., N, 
(~l) n 

P- 2 n(q) = —j=-Q n {q 2 - a - l]a,P,N), 

. (—1 )” 1 / (n + a + 1 )(n + a + (5 + l)( 2 n + 2 + a + f3) 

(ot - 1- 1) y (n + N + a + /3 + 2)(2n + a + /3 + 1) 

x qQn(q 2 -a-l;a + l,/3,N). 


These polynomials satisfy the discrete orthogonality relation 



(N - q 2 + a + /3 + 1\ 
\ AT — g 2 + a + 1 ) 


Pn(q)Pn'{q ) 


^"L n /2J (®’ Pi N)f3 n ,n' 


with 


S — | — \/ N -f Q! -)-1, — A ,r A n,..., — \ZoT-i-~"l, \/ o T^l 5 • • •, \/ A r -)- ot, y/ N -foil lj 

and 

h (n, R N\ = (~ 1 )"( ra + « + P + 1)7V+1 (/^ + 

(2n + a +/3 + l)(a + l) n (-Ai) n AM 

So P n (q) is a polynomial of degree n in the variable g, of different type (with different 
parameters when expressed as a Hahn polynomial) depending on whether n is even or n is odd. 
The support points of the discrete orthogonality are given by 

q = ±\/k + a + l, k = 0,..., N. 


In the same way, the dual orthogonality for the case Hahn II gives rise to 


Proposition 10. Let a > —1, j3 > —1, and consider the 2 N + 1 polynomials 


P- 2 n{q) = 


p2n+l{q) 


f—ll n 

—j=~Q n (q 2 -,a,P,N), n = 0,l,...,N, 

(— l) n 1 I (N — n)(n + a + l)(n + a + f3 + 1)(2 n + a + /3 + 2) 

\/2 (a + 1 ) A^ y ( 2 n + a + /? + 1 ) 

x qQ n {q 2 - l;« + l,/3,N - l), n = 0 , l,..., N — l. 


These polynomials satisfy the discrete orthogonality relation 

g ( q2 q2 a ) { N ~ ( [yy^ + Ko)Pn(q)Pn’(q) = hy n/ 2\ ( a i Pi N)Pn,ri 

with 


s= {-Vn,-Vn^i, ... ... ,Vn^i,Vn} 

and 

, / at, _ (-l) ra (n + a + ft + l)jy + i(/3 + l) n n! 

nl ,P ’ (2n + a + /3 + l)(a + l) n (-AT) n AT! ' 

The ideas described in the three propositions of this section should be clear. It would lead us 
too far to give also the explicit forms corresponding to the remaining cases. Let us just mention 
that also for these cases the support of the new polynomials coincides with the spectrum of the 
corresponding two-diagonal matrix M. 
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6 First application: eigenvalue test matrices 


In Sections 3 and 4 we have encountered a number of symmetric two-diagonal matrices M with 
explicit expressions for the eigenvectors and eigenvalues. In general, if one considers a two- 
diagonal matrix A of size (m + 2) x (m + 2), 


/ 0 b 0 0 \ 

c 0 0 h 

0 ci 0 ’' • 0 > 

bm 

\ 0 C m 0 J 


( 6 . 1 ) 


then it is clear that the characteristic polynomial depends on the products 6 jCj, i = 0 ,..., m, 
only, and not on bi and c, separately. So the same holds for the eigenvalues. Therefore, if all 
matrix elements bi and c* are positive, the eigenvalues of A or of the related symmetric matrix 


A' 


( 0 \/5oco 0 
Vboco 0 Vhci 
0 Vbici 0 




0 


V 




\/b r 


0 


are the same. The eigenvectors of A' are those of A after multiplication by a diagonal matrix 
(the diagonal matrix that is used in the similarity transformation from A to A'). 

For matrices of type (6.1), it is sufficient to denote them by their superdiagonal [b] = 
[bo,... ,b m ] and their subdiagonal [c] = [co,...,c m ]. So the Sylvester-Kac matrix from the 
introduction is denoted by 


[b] = [1, 2, ..., N], [c] = [IV, ..., 2,1], 

with eigenvalues given by ( 1 . 2 ). 

The importance of the Sylvester-Kac matrix as a test matrix for numerical eigenvalue routines 
has already been emphasized in the Introduction. In this context, it is also significant that the 
matrix itself has integer entries only (so there is no rounding error when represented on a digital 
computer), and that also the eigenvalues are integers. Of course, matrices with rational numbers 
as entries suffice as well, since one can always multiply the matrix by an appropriate integer 
factor. 

Let us now systematically consider the two-diagonal matrices encountered in the classification 
process of doubling Hahn or dual Hahn polynomials. For the matrix (3.18) of the dual Hahn I 
case, the corresponding non-symmetric form can be chosen as the two-diagonal matrix with 


[b] = [7 +1,1, 7 + 2, 2, ...,7 + AT, IV], 

[c] = [N,N + 5,N -1,N -1 + 6,...,1,5 + 1}. (6.2) 

The eigenvalues are determined by Proposition 2 and given by 0, ±y/k(k + 7 + 5 + 1), k = 
1 This is (up to a factor 2) the matrix (1.3) mentioned in the Introduction. As test 

matrix, the choice 7 + <5 + 1 = 0 (leaving one free parameter) is interesting as it gives rise to 
integer eigenvalues. In Proposition 2 there is the initial condition 7 > — l,h>— 1. Clearly, if 
one is only dealing with eigenvalues, the condition for (6.2) is just 7 + 5 + 2 > 0. And when 
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one substitutes <5 = —7 — 1 in ( 6 . 2 ), there is no condition at all for the one-parameter family of 
matrices of the form ( 6 . 2 ). 

For the dual Hahn II case, the matrix (3.22) is given in Proposition 3, and its non-symmetric 
form can be taken as 


[b] = [7 + AT, 1,7 + IV — 1 , 2 ,... ,7+ 1 ,IV], 

[c] = [N,S + l,N-l ,6 + 2,...,l,8 + N]. (6.3) 

The eigenvalues are given by 0, ±\/fc( 7 + 6 + 1 + 2 N — k), k = 1,... ,N. There is no simple 
substitution that reduces these eigenvalues to integers. 

For the dual Hahn III case, the matrix (2.15) is given in Proposition 4, and its simplest 
non-symmetric form is 

[b] = [7 +1,1,7+ 2, 2, . ..,7 +IV, IV, 7 +IV+ 1 ], 

[c] = [5 + N + 1, N,5 + N, N — 1,..., 5 + 2,1, <5 + 1], (6.4) 

The eigenvalues are given by (2.17), i.e., =t 7 / (7 -p/c-|- 1)(5 -p/c-p 1), k = 0,... , IV. Up to a fac¬ 
tor 2, this is the third matrix mentioned in the Introduction. The substitution 5 = 7 leads to 
a one-parameter family of two-diagonal matrices with square-free eigenvalues. And in particular 
when moreover 7 is integer, all matrix entries and all eigenvalues are integers. 

The two-diagonal matrices arising from the Hahn doubles or the Racah doubles can also be 
written in a square-free form of type (6.1). However, for these cases the entries in the two- 
diagonal matrices M are already quite involved (see, e.g., Propositions 6 , 7, 12 or 13), and we 
shall not discuss them further in this context. The three examples given here, (6.2)-(6.4), are 
already sufficiently interesting as extensions of the Sylvester-Kac matrix as potential eigenvalue 
test matrices. 


7 Further applications: related algebraic structures 
and finite oscillator models 


The original example of a (dual) Hahn double, described here in Section 2, was encountered in 
the context of a finite oscillator model [14]. In that context, there is also a related algebraic 
structure. In particular, the two-diagonal matrices M of the form (2.14) or (3.18) are interpreted 
as representation matrices of an algebra, which can be seen as a deformation of the Lie alge¬ 
bra su( 2 ). Once an algebraic formulation is clear, this structure can be used to model a finite 
oscillator. The close relationship comes from the fact that for the corresponding finite oscillator 
model the spectrum of the position operator coincides with the spectrum of the matrix M. 

Therefore, it is worthwhile to examine the algebraic structures behind the current matrices M. 
We shall do this explicitly for the three double dual Hahn cases. 

For the case dual Hahn I, we return to the form of the matrix M given in (3.18) or (3.19). 
For any positive integer N, let denote the lower-triangular tridiagonal ( 2N + 1) x (21V + 1) 
matrix given below, and let «/_ be its transpose 


/ 0 0 \ 
M 0 0 0 

0 Mi 0 0 

0 M 2 0 



(7.1) 


V 
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Let us also define the common diagonal matrix 


Jo = diag(—iV, —TV + 1,..., N), (7.2) 

and the “parity matrix” 

P = diag(l, —1,1, —1,...). (7.3) 

Then it is easy to check that these matrices satisfy the following relations (as usual, / denotes 
the identity matrix) 

P 2 = 1 , PJo = JoP, PJ± = ~J±P, 

[Jo, J±\ = ±J±, 

[J+, J-] = 2Jo + 2(7 + 5 + 1 ) J 0 P - (2N + 1)( 7 - 8 )P + (7 - 8)1. (7.4) 

Especially the last equation is interesting. From the algebraic point of view, it introduces some 
two-parameter deformation or extension of su(2). When 7 = 5 = —1/2, the equations coincide 
with the su(2) relations. Another important case is when 5 = —7 — 1, leaving a one-parameter 
extension of su( 2 ) without quadratic terms. 

For the case dual Hahn II, the corresponding expressions of J+, J_, Jo and P are the same 
as above in (7.1)—(7.3) , but with M^-values given by (3.22). As far as the algebraic relations are 
concerned, they are also given by (7.4) but with the last relation replaced by 

[J+, J_] = -2J 0 + 2(7 + 5 + 2N + 1) J 0 P + (2 N + 1 )( 7 - 5)P - (7 - 8)1. 

For the case dual Hahn III, the size of the matrices changes to (2A~ + 2) x (2 N + 2). For J + 
and J_ one can use (7.1), with M^-vallies given by (3.23). P has the same expression (7.3), but 
for Jo we need to take 

Jo = diag N — -, —N + N + -^ • 

With these expressions, the algebraic relations are given by (7.4) but with the last relation 
replaced by 


[J+, J_] = 2J 0 + 2(7 - 8 )J 0 P - (( 2 N + 2 )( 7 + <5 + 1) + (2 7 + 1)(2£ + 1 ))P 

+ (7 -8)1. (7.5) 

The structure of these algebras is related to the structure of the so-called algebra P, of the 
dual —1 Hahn polynomials, see [11, 31]. It is not hard to verify that the algebra P, determined 
by [11, equations (3.4)—(3.6)] or [11, equations (6.2)—(6.4)], can be cast in the form (7.4) (or vice 
versa). Indeed, starting from the form [11, equations (6.2)—(6.4)] coming from dual —1 Hahn 
polynomials, we can take 

Jo = Ki - J+ = K 2 + K 3 , J-=K 2 - K 3 , 

to get the same form as (7.4) 

P 2 = 1, PJo = JoP, PJ± = ~J±P , 

[J 0 , J±] = ±J±, [J+,J_] = 2Jo + 2 l /J 0 P+|p+^I, (7.6) 

where v, a , p depend on the parameters of the dual —1 Hahn polynomials a, j3, N through [11, 
equations (3.4)-(3.6)]. In our case, the algebraic relations are the same, but the dependence of 
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the “structure constants” in (7.6) on the parameters 7 , 6 , N of the dual Hahn polynomials is 
different. 

As far as we can see, the doubling of dual Hahn polynomials as classified in this paper 
gives a set of polynomials that is similar but in general not the same as a set of dual — 1 
Hahn polynomials [31] (except for specific values of parameters, e.g., <5 = —7 — 1 does coincide 
with a specific dual —1 Hahn polynomial). For general parameters, the support of the weight 
function is different, the recurrence relations (or difference relations) are different, and the 
hypergeometric series expression is different. 

The algebraic structures obtained here (or special cases thereof) can be of interest for the 
construction of finite oscillator models [1, 2, 3, 14]. Two familiar finite oscillator models fall 
within this framework: the model discussed in [14] corresponds to (7.5) with (5 = 7 , and the 
one analysed in [15] to (7.4) with <5 = 7 . Observe that there are some other interesting special 
values. For example, the case (7.4) with 6 = —7 — 1 gives rise to an interesting algebra, and in 
particular also to a very simple spectrum (3.20). We intend to study the finite oscillator that 
is modeled by this case, and study in particular the corresponding finite Fourier transform; but 
this will be the topic of a separate paper. 

8 Conclusion 

We have classified all pairs of recurrence relations for two types of dual Hahn polynomials (i.e., 
dual Hahn polynomials with different parameters), and refer to these as dual Hahn doubles. 
The analysis is quite straightforward, and the result is given in Theorem 1, yielding three cases. 
For each case, we have given the corresponding symmetric two-diagonal matrix M, its matrix 
of orthonormal eigenvectors U and its eigenvalues in explicit form. The same classification has 
been obtained for Hahn polynomials and Racah polynomials. 

The orthogonality of the matrix U gives rise to new sets of orthogonal polynomials. These 
sets could in principle also be obtained from, for example, a set of dual Hahn polynomials and 
a certain Christoffel transform. In our approach, the possible cases where such a transform gives 
rise to a polynomial of the same type follow naturally, and also the explicit polynomials and 
their orthogonality relations arise automatically. 

As an interesting secondary outcome, we obtain nice one-parameter and two-parameter ex¬ 
tensions of the Sylvester-Kac matrix with explicit eigenvalue expressions. Such matrices can be 
of interest for testing numerical eigenvalue routines. 

The first example of a (dual) Hahn double appeared in a finite oscillator model [14]. For this 
model, the Hahn polynomials (or their duals) describe the discrete position wavefunction of the 
oscillator, and the two-diagonal matrix M lies behind an underlying algebraic structure. Here, 
we have examined the algebraic relations corresponding to the three dual Hahn cases. It is clear 
that the analysis of finite oscillators for some of these cases is worth pursuing. 


A Appendix: doubling Racah polynomials 


The technique presented in Sections 3 and 4 is applied here for Racah polynomials. 

Racah polynomials R n (X(x); a, /?, 7 , <5) of degree n (n = 0,1,..., N ) in the variable X(x) = 
x(x + 7 + 6 + 1) are defined by [13, 19, 23] 


R n (A(x);a,/3,7, <5) = 4 F 3 


—n, n + a + /3 + 1, —x, x + 7 + <5 + 1 


a + l,/3 + <5 + l, 7 Tl 
where one of the denominator parameters should be — N: 


a + 1 = — N or 


f3 + 6 + 1 = -N 


or 


7 + 1 = —N. 


(A.l) 
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For the (discrete) orthogonality relation (depending on the choice of which parameter relates 
to — N ) we refer to [19, equation (9.2.2)] or [25, Section 18.25] 

Racah polynomials satisfy a recurrence relation of the form (2.4) with 


y n {x) = i?„(A(x);a,/ 3 , 7 , 5 ), A(x) = A(x) = x(x + 7 + 5 + 1), 
. . (n + ot ++)(n + 01 + f3 + \)(n +'y + T)(n + j3 + 5 + \) 

(2 n + a + (3 + l)( 2 ?r + a + (3 + 2 ) 

r( _ n ( n + a + P ~'y)( n + a ~ ^)( n + P) 

(2n + a + f3)(2n + a + (3 + 1) 


(A.2) 


We have applied the technique described in (3.3)-(3.15), with y n = R n (A(x); a, j3, 7 , 5) and 
y n = R n (A(x); a, f3, 7 , 5). The analysis is again straightforward but tedious, and the final out¬ 
come is 


Theorem 11. The only way to combine two sets of Racah polynomials such that they satisfy 
difference relations of the form (3.3), (3.4) is one of the four cases: 

Racah I, R n (x ) = R n ( A(x); a, (3, 7 , <5) and R n (x) = R n ( A(x); a, (3 + 1 ,7 + 1, <5 — 1): 

(n + /3 + 6 + 1 )(n + ol + 1) , , ( n — 6 + ol + l)(n + f3 + 1) , . 

-R n+ i(x) -——■- ' - Rn{x) 


(2 n + a + (3 + 2) 

(x + 5)(x + 7 + I) 


(2n + ol + (3 + 2) 


7+1 

(n + a + (3 + 2 ) (n + 7 + 2 ) 
(2n + a + ft + 3) 

= (7+ l)Rn+l(x). 


Rnfx), 


k _ (n+lHn-T + a + l + l) ^ 
+u ; (2n + a + /3 + 3) v ; 


Racah II, R n (x) = R n (X(x); a, (3, 7 , (5) andR n (x) = R n (A(x);a,/3 + l, 7 , 5 ): 

(n + 7 +l)(n + a + 1 ) (n - 7 + a +/3 + l)(n +/? + 1 ) 

-R n+ i(x)-—— - -Rn(^) 


( 2 n + a + /3 + 2 ) 

(x + (3 + <5 + 1) (x + 7 - /3) 


( 2 n + 01 + /3 + 2 ) 


(3 + <5 + 1 

(n + f3 + d + 2 ) (?r + o + /3 + 2 ) 
(2n + a + (3 + 3) 

= (/3 + 5 + l)R n +i(x). 


Rn ( X ), 


(„+!)(„-S + a+D ^ 

+ v ’ (2n + a + /3 + 3) v 7 


Racah III, R n (x) = i? ra (A(x); a, (3, 7 , <5) and R n (x) = R n ( A(x — 1); a + 1, /3 ,7 + 1, <5 + 1): 

1 „ 1 


Rn+l(®) 


n / \ _ x(x + 7 + (5 + 1) ~ 

Rn\X) — -i\^n\X)i 


(2n + a + f3 + 2 ) ( 2 n + a + (3 + 2 ) (7 + l)(/3 + (5 + l)(a + 1) 

(71 + 7 + 2 )(n + /3 + 5 + 2)(n + a + 2)(n + a + (3 + 2) 


/r) , . Q 1 Q\ Rn-fl(+) 

(2n -\- ql -\- p -\- 3) 

(n + l)(n - 7 + a + f3 + l)(n - <5 + a + l)(n + /3 + 1) 
(2n + ol + f3 + 3) 

= (7 + 1 )(^ + <5 + l)(a + l)R n +i(x). 


R n (x) 


Racah IV, R n (x) = R n ( A(x); a, f3, 7, 5) and R n (x) = i?„(A(x); a + 1 ,(3, 7, 5 ): 

(n + 7 + l)(n + f3 + S + 1) . . (n - 7 + a + f3 + l)(x — S + a + 1) . . 

■Rn+lV.X) — | | — ■ —r R n [x) 


(2n + a. + (3 + 2) 


( 2 n + o + f3 + 2 ) 
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_ (x + 7 + 5 - a)(x + a + 1) £, ^ 

— 7 777 Rn\ x )i 

(a + 1) 

[n T ol + 2) (n + a + P -P 2) ~ , . (u + 1) (??, + f3 + 1) ~ . . 

/O I , /i q\ -TL^-)-l (XJ . i /3 i o\ t C \ x ) 

[f2n -\- ot ~\~ p ~\~ 3) (2u 1“ o: h p 4" 3) 


(a + l)i 4 + i(x). 


Note that after interchanging n and x, and a 0 7 and /?•<->• 5, the relations in Racah III 
coincide with the known forward and backward shift operator relations [19, equations (9.2.6) 
and (9.2.8)]. The relations in Racah I were already found in [16, equations (5) and ( 6 )]. 

In the context of Section 5 it is worth noting that the above relations also correspond to 
Christoffel-Genonimus transforms. Taking R n (x) = R n (\(x)',a,/3,j,8) in the relations (5.3)- 
(5.5), with data given by (A.2), the above cases Racah I, II, III, IV correspond respectively 
to the choices 12 = —5, 12 = (5 — 7 , v = 0 and v = —a — 1. 

For each of the four cases, one can translate the set of difference relations to a matrix identity 
of the form MU = UD. In fact, for each of the four cases, there are three subcases depending on 
the choice of —N in (A.l). We shall not give all of these cases: they should be easy to construct 
for the reader who needs one. Let us just give an example or two. 

Consider the case Racah I with a + 1 = —N. It is convenient to perform the shift <5 —> 
5 + 1 in the two difference relations of Theorem 11. The orthogonal matrix U is of order 
(2 N + 2 ) x (2 N + 2 ), with elements 


U'2n,N—x — U2n,N+x+\ — 


(- 1 )’ 

V2 


■-Rn(A(x); a, j3, 7 ,5 + 1 ), 


U2n+l,N—x — U2n+l,N+x+l — 


V2 


-Rn(X(x)-,a,/3 + 1,7 + 1,<5), 


(A.3) 


where R n is the notation for a normalized Racah polynomial. Then, one has 

Proposition 12. Suppose that 7 ,5 > —1 and /3 > N + 'y or /3 < —N — 5—1. Let M be 
a tridiagonal (2 N + 2) x (2 N + 2 )-matrix of the form (2.14) with 



(N - (3 — k)( 7 + 1 + k)(N + 5 + 1 - k)(k + /3 + 1) 
(N - p - 2k)(2k - N + l + P) : 


(7 + N - p - k){k + 1)(N - k)(k + P + <5 + 2) 
(N - p - 2k - 2){2k - N + 1 + p) 


and U the orthogonal matrix determined in (A.3). Then the columns of U are the eigenvectors 
of M, i.e., MU = UD, where D is a diagonal matrix containing the eigenvalues of M 


D = diag(—ejv,..., -ei, -e 0 , £ 0 , 6 !,..., e N ), 

£fc = \J(k + 7 + l)(fc + <5 + 1), k = 0,1,..., N. 

As a second example, consider the case Racah III with a + 1 = —N. The orthogonal 
matrix U is now of order ( 2 N + 1 ) x ( 2 N + 1 ), with elements 

(-l) n ~ 

U‘2n,N—x — U 2 n,N+x — PiTi ^i ^ AT, 

(-l) n - 

U2n+i,N-x-i = ~ U2n+i,N+x+i =- R n ( A(x );a + 1,/3,7 + 1 ,5 + 1), 

n = 0,..., N — 1, 

U 2n ,N = (-ir^n(X(0);a,P, 7 ,5), U 2n+ i,N = 0. (A.4) 


Then, one has 
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Proposition 13. Suppose that 7 ,5 > —1 and (3 > N + 7 or (3 < —N — 5. Let M be a tridiagonal 
(2 N + 1) X (2N + 1)-matrix of the form (3.18) with 



(k + 7 + 1)(—iV + (3 + k)(N — k)(k + f3 + 6 + 1) 
{N - p - 2k)(N - p - 2k - 1) 


(7 + N — P — k)(k + l)(k + (3 + 1)(A: — 5 — N) 
(N — P — 2k — 2)(N -p-2k-l) : 


and U the orthogonal matrix determined in (A.4). Then the columns of U are the eigenvectors 
of M, i.e., MU = UD, where D is a diagonal matrix containing the eigenvalues of M 

D = diag(-ejv,...,-ei,0,ei,...,ejv), e k = \/k(k + 7 + 5 + 1), k = l,...,N. 
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